Violation of Entropic Leggett-Garg Inequality in Nuclear Spins 
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We report an experimental study of recently formulated entropic Leggett-Garg inequality (ELGI) 
by Usha Devi et al. (arXiv: 1208.4491v2 (2012)). This inequality places a bound on the statistical 
measurement outcomes of dynamical observables describing a macrorealistic system. Such a bound is 
not necessarily obeyed by quantum systems, and therefore provides an important way to distinguish 
quantumness from classical behavior. Here we study ELGI using a two-qubit nuclear magnetic 
resonance system. To perform the noninvasive measurements required for the ELGI study, we 
prepare the system qubit in a maximally mixed state as well as use the 'ideal negative result 
measurement' procedure with the help of an ancilla qubit. The experimental results show a clear 
violation of ELGI by over four standard deviations. These results agree with the predictions of 
quantum theory. The violation of ELGI is attributed to the fact that certain joint probabilities 
are not legitimate in the quantum scenario, in the sense they do not reproduce all the marginal 
probabilities. Using a three-qubit system, we experimentally demonstrate that three-time joint 
probabilities do not reproduce certain two-time marginal probabilities. 
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Introduction. — The behavior of quantum systems is of- 
ten incomprehensible by classical notions, the best exam- 
ples being nonlocality [1, 2] and contextuality [3]. Quan- 
tum systems are nonlocal since they violate Bell's in- 
equality, which assumes that local operations on one of 
the two space-like separated objects can not disturb the 
measurement outcomes of the other [4]. The quantum 
systems are also contextual in the sense that a measure- 
ment outcome depends not only on the system and the 
property being measured, but also on the context of the 
measurement, i.e., on the set of other compatible prop- 
erties which are being measured along with. 

Another notion imposed on classical objects is macro- 
realism, which is based on two criteria: (i) the object 
remains in one or the other of many possible states at 
all times, and (ii) the measurements are noninvasive, i.e., 
they reveal the state of the object without disturbing 
the object or its future dynamics. Quantum systems 
are incompatible with these criteria and therefore vio- 
late bounds on correlations derived from them. For in- 
stance, Leggett-Garg inequality (LGI) sets up macrore- 
alistic bounds on linear combinations of two-time corre- 
lations of a dichotomic observable belonging to a single 
dynamical system [5]. In this sense, LGI is considered as 
a temporal analogue of Bell's inequality. Quantum sys- 
tems do not comply with LGI, and therefore provide an 
important way to distinguish the quantum behavior from 
macrorcalism. Violations of LGI by quantum systems 
have been investigated and demonstrated experimentally 
in various systems [6-9]. 

For understanding the quantum behavior it is impor- 
tant to investigate it through different approaches, par- 
ticularly from an information theoretical point of view. 
For example, an entropic formulation for Bell's inequality 



has been given by Braunstcin and Caves [10], and more 
recently that for contextuality has been given indepen- 
dently by Rafael and Fritz [11] and Kurzyhski et.al. [12]. 
Recently, an entropic formulation of LGI has also been 
introduced by Usha Devi et al. [13], in terms of classical 
Shannon entropies associated with classical correlations. 
Such entropies obey certain constraints, which when vio- 
lated would imply non-existence of legitimate joint prob- 
abilities (JP) for all the measurement outcomes. 

Here we report an experimental demonstration of vi- 
olation of entropic LGI (ELGI) in an ensemble of spin 
1/2 nuclei using nuclear magnetic resonance (NMR) tech- 
niques. Although NMR experiments are carried out at a 
high temperature limit, the nuclear spins have long co- 
herence times, and their unitary evolutions can be con- 
trolled in a precise way. The large parallel computa- 
tions carried out in an NMR spin ensemble assists in effi- 
ciently extracting the single-event probability (SEP) and 
JP. The simplest ELGI study involves three sets of two- 
time joint measurements of a dynamic observable belong- 
ing to a 'system' qubit at time instants (£1,^2), (£2,^3), 
and (ti , £3). The first measurement in each case must 
be 'noninvasive' in the sense, it should not influence the 
outcome of the second measurement. These noninvasive 
measurements (NIM) can be performed with the help of 
an ancilla qubit. 

Further, it has been argued in [13] that the violation of 
ELGI arises essentially due to the fact that the JP do not 
originate from a legitimate grand probability (of which 
the JP are the marginals). Here we describe extracting 
three-time JP using a three-qubit system, and demon- 
strate experimentally that it can not reproduce all the 
marginal probabilities (MP) substantiating this feature. 

In the following we briefly revisit the theory of ELGI 
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[13] and then we describe the circuits for the measure- 
ment of SEP and JP. Later we detail its experimental 
study using a two-qubit NMR system. Then we describe 
the study of three-time JP using a three-qubit NMR sys- 
tem. 

Theory. — Consider a dynamical observable Q{tk) = 
Qk measured at different time instances t k ■ Let the mea- 
surement outcomes be qk with probabilities P(q k )- In 
classical information theory, the amount of information 
stored in the random variable Qk is given by the Shannon 
entropy [14], 



H(Qk) = -^2 P ( ( lk)\og 2 P(qk). 



(1) 



The conditional information stored in Qk+i at time tk+i , 
assuming that the observable Qk has an outcome qk, is 

H(Qk+i\Qk = ?fc) = - X] P (<lk+i\qk)\og 2 P(qk + i\qk), 

1k + l 

where P(qk+i \qk) is the conditional probability. Then the 
mean conditional entropy is given by, 

H(Q k +i\Qk) = -Y, P MH(Qk+i\Qk^qk)- (2) 



Using Bayes' theorem, P(q k+ i\qk)P(qk) = P(qk+l,qk), 
the mean conditional entropy becomes 

H(Qk+i\Qk) = H(Qk, Qk+i) - H(Q k ), (3) 
where the joint Shannon entropy is given by 

H(Q k ,Qk+i) = - P(qk+l,qk)log 2 P(q k+ i,q k )(4) 

IkAk + l 

These Shannon entropies always follow the inequality [10] 

H(Q k+l \Q k ) < H(Q k+ i) < H(Q k ,Q k+ i). (5) 

The left side of the equation implies that removing a con- 
straint never decreases the entropy, and the right side im- 
plies information stored in two variables is always greater 
than or equal to that in one [13] . Suppose that three mea- 
surements Qk, Qk+i, and Qk+ m , are performed at time 
instants t k < t k +i < tk+ m - Then, from equations (3) and 
(5), the following inequality can be obtained: 

H{Q k+m \Qk) < H{Qk +m \Qk+i) + H(Q k+ i\Qk). (6) 

For n measurements Q\, Q2, ■ ■ ■ , Q n , at time instants 
<i < t 2 < ■ ■ ■ < t„, the above inequality can be gen- 
eralized to [13] 



J2 H{Qk\Qk-i) - H(Qn|Qi) > 0. 



(7) 



k=2 



This inequality must be followed by all macro-realistic 
objects, since its satisfaction means the existence of le- 



gitimate JP distribution, which can yield all MP [12]. 

Usha Devi et al. [13] have shown theoretically that 
the above inequality is violated by a quantum spin-s 
system, prepared in a completely mixed initial state, 
Pin = l/(2s + 1). Consider the z-componcnt of the 
spin evolving under the Hamiltonian % = —luS x as 
our dynamical observable, i.e. Qt = UtS z U}, where 
Ut = e~ lHt , and S x and S z are the components of spin- 
angular momentum. Let n-measurements occur at regu- 
lar time instants At, 2At, ■ ■ ■ , nAt. Ideally in this case, 
the conditional entropies H (Qk\Qk~i) between succes- 
sive measurements are all equal, and can be denoted as 
H[9/(n — 1)], where 9/(n — 1) = a; At is the rotation 
caused by the Hamiltonian in the interval At. Similarly 
we can denote H(Q n \Qi) as H[0\. The lhs of inequality 
(7) scaled in units of log 2 (2,s + 1) is termed as the in- 
formation deficit T>. For n-cquidistant measurements, it 
can be written as [13] 



V n {9) 



(n-l)H[9/(n-l)]-H[6] 
log 2 (2s + l) 



> 0. 



(8) 



Measurement of Probabilities. — Consider a spin-1/2 
particle as the system qubit. Using the eigenvectors 
{|0), |1)} of S z , as the computational basis, the projec- 
tion operators at time t — arc {II Q = \a) (a|} Q =o,i- 
For the dynamical observable, the measurement basis 
is rotating under the unitary Ut = e iwSxt , such that 
= UJlaU}. However, it is convenient to perform 
the actual measurements in the time-independent com- 
putational basis. Since for an instantaneous state p(t), 

n^p(t)n^ = U t U a (u}p(t)Utj U a u}, measuring in {n^} 

basis is equivalent to back-evolving the state by If} , mea- 
suring in computational basis, and lastly forward evolv- 
ing by U t . This latter evolution can be omitted if one is 
interested only in the probabilities and not in the post 
measurement state of the system. For example, in case of 
multiple-time measurements, the forward evolution can 
be omitted after the final measurement. 

The method for extracting SEP and JP involves the 
quantum circuits shown in Fig. 1. To measure SEP 
P(qi) of system qubit in a general state ps, it is evolved 



byUj 



and the probabilities P{qi) are obtained 



using diagonal tomography (Fig. 1(a)). Here a further 
forward evolution by Ui is not necessary as described ear- 
lier. 

To measure JP P(qi,qj), we utilize an ancilla qubit 
initialized in the state |0)(0| (Fig. 1(b)). After back 
evolution to computational basis, the CNOT gate en- 
codes the probabilities of the system-qubit P(qi) on 
to the ancilla-qubit [15]. After a further evolution by 
UiUj = e ~ lw5 '^( t J- t i) j a diagonal tomography of the two 
qubit system yields P(qi,qj) [15]. 

A similar scheme, shown in Fig. 1(c), is employed for 
extracting three-time JP. These circuits can be general- 
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FIG. 1. Circuits for measuring SEP (a), two-time JP (b), and 
three-time JP (c). The grouped gates represent measurement 
in {IIo, 111} basis. In (b) the operation C can be either CNOT 
or anti-CNOT gate as described in the text. 
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FIG. 2. The lines indicate theoretical JP and the symbols in- 
dicate the mean experimental probabilities obtained by INRM 
procedure. 



ized for higher order JP or for spin-s > 1/2 systems, 
using appropriate ancilla register. 

In the earlier LGI experiments, NIM have been per- 
formed by either (i) a weak measurement which causes 
minimum disturbance to the quantum state [7] or (ii) 
initializing the system qubit in a maximally mixed state 
so that the system density matrix remains unchanged by 
the measurements [8]. Recently however, it was noted by 
Knee et al. that a sceptical macrorealist is not convinced 
by either of the above methods [16]. Instead, they had 
proposed a 'ideal negative result measurement' (INRM) 
procedure that is more convincingly noninvasive [9] . The 
idea is as follows. The CNOT gate is able to flip the 
ancilla qubit only if the system qubit is in state |1), and 
does nothing if the system qubit is in state |0). There- 
fore after the CNOT gate, if we measure the probability 
of unflipped ancilla, this corresponds to an 'interaction- 
free' or NIM of P(q = 0). Similarly, we can implement 
an anti-CNOT gate, which flips the ancilla only if the 
system qubit is in state |0), and docs nothing otherwise, 
such that the probability of unflipped qubit now gives 
P(q = 1). Note that in both the cases, the probabilities 
wherein the system interacted with the ancilla, resulting 
in its flip, arc discarded. The final measurement need 
not be NIM since we are not concerned about any fur- 
ther evolution. 

In our experiments we combine the two methods, i.e., 
(i) first we prepare the system in a maximally mixed state 
i.e., ps = 1/2, and (ii) we perform INRM. In this case, 
F(0i) = P(U) = 1/2, and JP are 

P(0,,0j) = |cos(%/2)| 2 /2 = P{U,l ), and, 
P(0 i ,l J ) = |sin(%/2)| 2 /2 = P(li,0 i ) J (9) 

where 9ij = ui(tj — U) [13]. 

The only single event entropy needed for the ELGI 
test is H(Qi), since H(Q t ) is constant for the maxi- 
mally mixed system state. Further, since H(Qi,Q2) — 



H(Q2, Qz) in the case of uniform time intervals, only two 
joint entropies H(Qi,Q 2 ) and H(Qi,Q 3 ) are needed to 
be measured for evaluating D3. In the following we de- 
scribe the experimental implementation of these circuits 
for the three- measurement LGI test. 

Experiment. — We have used 13 CHCi3 (dissolved in 
CDCI3) as the two qubit system and treat its 13 C and 
*H nuclear spins as the system and the ancilla qubits re- 
spectively. The resonance offset of 1 C was set to 100 Hz 
(provides a dynamic 'observable') and that of 1 H to Hz 
(on resonant). The two spins have an indirect spin-spin 
coupling constant J = 209.2 Hz [15]. All the experiments 
were carried out at an ambient temperature of 300 K on 
a 500 MHz Bruker UltraShield NMR spectrometer. 

The evolution propagator Uj^Ui = g-^^w-U) j s re _ 
alized by the cascade HE/<jlrI, where H is the Hadamard 
gate, and the delay propagator Ud = e~ %SzU! ^ ti ~ t ^ > cor- 
responds to the z-precession of the system qubit at 
uj = 27rl00 rad/s resonance off-set. The J-evolution dur- 
ing this delay is refocused by a 7r pulse on the ancilla 
qubit. The CNOT, H, as well as the ir pulses are realized 
by numerically optimized amplitude and phase modu- 
lated RF pulses, and are robust against the RF inho- 
mogencity with a average Hilbcrt-Schmidt fidelity better 
than 0.998 [17-19]. The final measurement of probabili- 
ties are carried out by diagonal tomography. It involved 
dephasing all the coherences using a strong pulsed field 
gradient followed by a 7r/30 detection pulse. The inten- 
sities of the resulting spectral lines yielded a traceless di- 
agonal density matrix, which was normalized and added 
with identity matrix to extract the probabilities. As de- 
scribed in Fig. 1(b), two sets of experiments were per- 
formed, one with CNOT and the other with anti-CNOT. 
We extracted P(0,q) (q = {0,1}) from the CNOT set, 
and P(l,q) from the anti-CNOT set. The probabilities 
thus obtained by INRM procedure are plotted in Fig. 2. 
These sets of experiments also allow us to compare the 
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FIG. 3. Information deficit T>z versus 8 obtained with 
CNOT (a,b), anti-CNOT (c,d), and INRM procedure (e,f). 
The boxed areas in the left plots (a,c,e) are magnified in the 
right plots (b,d,f) respectively. The mean experimental T>z 
(in bits) values are shown as symbols. The curves indicate 
theoretical T>z (in bits). The horizontal line T>i = indicates 
the lower bound of the macrorealism territory. 



i.e., P'(q 1 ,q 2 ) = P{qi,q 2 ) and P'{q 2 ,qa) = P{q2,qs)- 
However, P{q\, q 2 , q 3 ) can not reproduce P{qi,q3), i.e., 
P'{<li,q.s) ^ P{Qi,Q3), m general [13]. 

The above concept can be investigated experimentally 
by measuring the three-time JP, as described in Fig. 1 (c) . 
Since this experiment requires measurements at three 
time instants, we need two ancilla qubits along with the 
system qubit. We use the three 19 F nuclear spins (spin- 
1/2) of trifiuoroiodoethylene dissolved in acetone-D6 as 
the three-qubit system [15]. The chemical shifts differ- 
ences are ^2 — ^1 = 11860.8 Hz, v 2 — v 3 = 17379.1 Hz, 
and the scalar coupling constants are: J12 = 69.9 Hz, 
J13 = —128.3 Hz, and J23 = 47.4 Hz. Here the first spin 
(Fi) is used as the system qubit and the others (F 2 and 
-F3) are chosen as the ancilla qubits. Initialization in- 
volved preparing the state, ±=n + e {±l s ® |00)(00| j4 } 
where e ~ 10~ 5 is the purity factor [15, 20]. The ex- 
perimental three-time JP P{q\, q 2 , 93) obtained using the 
circuit Fig. 1(c) are shown in Fig. 4(a). Two-time JP 
P(qi,qj) were also obtained using a similar circuit (Fig. 
1(c) without Ujlll). Here JP are completely stored in the 
ancilla qubits and were obtained by tracing out the sys- 
tem qubit. The results P{qi,q 2 ), P(q 2 ,q 3 ), and P(qi,q 3 ) 
are shown in Fig. 4 (b-d). In each plots, we have over- 



results from (i) only CNOT, (ii) only anti-CNOT, and 
(iii) INRM procedures. The joint entropies were calcu- 
lated in each case using the experimental probabilities 
and the information deficit (in bits) was calculated us- 
ing the expression V 3 = 2H(Q 2 \Qi) - H{Q 3 \Q 1 ). The 
theoretical and experimental values of T> 3 for various ro- 
tation angles 9 are shown in Fig. 3. We find a gen- 
eral agreement between the mean experimental T> 3 val- 
ues with that of the quantum theory. The error bars 
indicate the standard deviations obtained by a series of 
independent measurements. According to quantum the- 
ory, a maximum violation of T> 3 = —0.134 should occur 
at 6 = 7t/4. The experimental values of D3(7r/4) are 
-0.141 ± 0.005, -0.136 ± 0.002, and -0.114 ± 0.027 for 
the CNOT, anti-CNOT, and INRM cases respectively. 
Thus in all the cases, we found a clear violation of ELGI. 

Three-time JP. — It has been argued that the viola- 
tion of ELGI in quantum systems is essentially because 
certain JP can not reproduce all MP [13]. For example, 
in the T> 3 experiment described earlier, the two-time JP 
P(qi,q 2 ), P(q 2 ,q 3 ), and P{q\,q 3 ) can be obtained from 
cqns. (9). From the three-time JP, it is possible to gen- 
erate MP: 

P' {li,Q.2) = Y,q 3 P{Qi,Q2,Qa), 
P'{Q3,Qs) = T, qi P{Qi,Q2,Q3), and 
P'{quq3) = E g2 P{Qi,q2,q 3 ). (10) 

Now P{qi,q 2 ,q 3 ) can reproduce P(qi,q 2 ) and P(q 2 ,q 3 ), 
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laycd MP P'(q 1 ,q 2 ), P'(q 2 ,q 3 ), and P'(qi,q 3 ) generated 
from the experimental three-time JP. As expected, MP 
P'io.XiQ.2) and P'(q 2 ,q 3 ) match well with the JP P(qi,q 2 ) 
and P(q 2 ,qz) respectively, while the MP P'{qi,qz) does 
not match with JP P(qi, q 3 ). These results confirm that 
the three-time JP P{qi,q 2 ,q 3 ) is not legitimate in the 
quantum case. It is interesting to note that even for 
those values of 9 for which 2?3 is positive, the three-time 
JP is illegitimate. Therefore, while the violation of ELGI 
indicates the quantumness of the system, its satisfaction 
docs not rule out the quantumness. 

Conclusions. — We described an experimental study 
of the entropic Leggett-Garg inequality in nuclear spins 
using NMR techniques. We employed the recently de- 
scribed 'ideal negative result measurement' procedure 
to noninvasively extract joint probabilities. Our re- 
sults indicate the macrorealistic bound being violated 
by over four standard deviations, confirming the non- 
macrorealistic nature of the spin-1/2 particles. Quantum 
systems do not have legitimate joint probability distribu- 
tion, which results in the violation of bounds set-up for 
macrorealistic systems. We have experimentally mea- 
sured the three-time joint probabilities and confirmed 
that the two-time joint probabilities are not reproduced 
as marginals. 

One distinct feature of the entropic LGI is that, the 
dichotomic nature of observables assumed in the original 
formulation of LGI can be relaxed, thus allowing one to 
study the quantum behavior of higher dimensional sys- 
tems such as spin > 1/2 systems. This could be an in- 
teresting topic for future experimental investigations. 
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SUPPLEMENTARY INFORMATION FOR "VIOLATION OF ENTROPIC LEGGETT-GARG 

INEQUALITY IN NUCLEAR SPINS" 



Encoding the probability using a CNOT gate: Consider a system qubit initially prepared in a general state 
ps and an ancilla qubit prepared in the state |0)(0|. The state of the system qubit can be written as P(0j)|0)(0| + 
P(lj)|l)(l| +a|l)(0| +aT|0)(l|, where a is a complex probability amplitude. The CNOT gate encodes the probability 
of the outcomes in the diagonal elements of ancilla qubit since, 

(P(0i)\0){0\ s + P(li)|l)(l| s + a|l)(0| s + a t |0)(l| s ) ® |0)<0|^ C ™ T |0)(0| S ® P(0 4 )|0)<0|^ + |1)(1| 5 ® P(1 4 )|1)(1| A 

+ |l)(0| s ®o|l)(0| A + |0)<l| s ®ot|0)<l| A . 

The probabilities P(0j) and P(lj) can now be retrieved by tracing over the system qubit and reading the diagonal 
elements of the ancilla state. 



The qubit systems: Fig. 5 shows the molecular structures and the Hamiltonian parameters of 13 CHCl3 (Fig. 
5(a,b)) and trifluoroiodoethylene (Fig. 5(c,d)). In the case of 13 CHCl3, spin-lattice (Tl) and spin-spin (T2) relaxation 
time constants for the 1 H spin are, respectively, 4.1 s and 4.0 s. The corresponding time constants for 13 C are 5.5 s 
and 0.8 s. In the case of trifluoroiodoethylene, the effective 19 F transverse relaxation time constants (T| ) were about 
0.8 s and their longitudinal relaxation time constants were all longer than 6.3 s. 
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FIG. 5. The molecular structures of chloroform (a) and trifluoroiodoethylene (c) and the corresponding tables (b and d) of 
relative resonance frequencies (diagonal elements) and the J-coupling constants. The pulse sequence for initializing trifluoroio- 
doethylene is shown in (e). In (e) the open pulses are tt pulses and the delay r = 1/(4 J23). 



Initialization in two-qubit ( 13 CHCls) system: The initialization involved preparing the maximally mixed state 
ps = 1/2 on the system qubit ( 13 C). This is achieved by a ir/2 pulse on 13 C followed by a strong Pulsed Field 
Gradient (PFG). 

Initialization in three-qubit (trifluoroiodoethylene) system: The equilibrium deviation density matrix evolves 
under the PPS sequence (Fig. 5(e)) as follows 

c . c , o W2)i,(V3)3x,PFG 2 1 (7T/4W 1 1 1 l/(2./ ?3 ) . /- , 

Jlz + J2z + J3z > ^z ' 2 3z 7/2 \/2 V 2 -^^2z + V -ZJ2x03z + 2^3z 

I (V4)- 2y ,PFG 

The above deviation density matrix is equivalent to the traceless part of -^pl + e { 5I5 ® |00)(00| j4 } where 
e ~ 10 -5 is the purity factor [20]. 

Diagonal tomography: 

The diagonal tomography was carried out using a strong PFG to dephase out all the residual off-diagonal elements 
and using a 6° non-selective linear readout pulse. The resulting intensities of the spectral lines constrain the diagonal 
elements (di) of the traceless deviation density matrix. The experimental deviation density matrix is normalized 
w.r.t. the theoretical traceless density matrix such that they both have the same root mean square value v'S?: > 
and trace is introduced by adding identity matrix to the normalized deviation matrix. The resulting diagonal den- 
sity matrix yields the probabilities. Estimation of random errors were carried out by several repetitions at each 9 value. 



